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Abstract 

We show that the hkehhood ratio of Higgs search experiments is a form to 
report the experimental results suitable to be combined with the information 
from precision measurements to obtain a joint constraint on the Higgs mass. 
We update our previous combined analysis using the new results on direct 
searches and recent precision measurements, including also the Z° leptonic 
partial width result. The method is also improved to take into account small 
non linearity effects in the theoretical formulae. We find an expected value 
for the Higgs mass around 160-170 GeV with an expectation uncertainty, 
quantified by the standard deviation of the distribution, of about 50-60 GeV. 
The 95% probability upper limit comes out to be around 260-290 GeV. 



*Contributed paper to the Workshop on "Confidence Limit", CERN, Geneva, Switzer- 
land, Jan. 17-18 2000. 



1 Introduction 



The Higgs boson is still the missing particle of the Standard Model (SM) 
picture and a considerable effort has been devoted to search for evidence 
of it. Unfortunately, till now all direct search experiments have been un- 
successful. However, the impressive amount of data collected at LEP, SLC, 
and the Tevatron allows to probe the quantum structure of the SM, thereby 
providing indirect information about the Higgs mass. While the negative 
outcome of the Higgs searches at LEP is usually reported as a combined 95% 
Confidence Level (C.L.) lower bound, the virtual Higgs effects are analyzed 
through a fit to the various precision observables that allows a 95% C.L. 
upper bound to be derived. In Ref.[||] we proposed a method to combine 
the information on the Higgs boson coming from direct searches (that we in- 
dicate generically as dir.) with that obtained from precision measurements 
(ind.), in order to derive a probability density function (p.d.f.) for its mass 

/(mH|"data","SM") = f{mH\ dir. k ind.) 

conditioned by both kind of experimental results under the assumption of 
validity of the SM. The heart of our method is the use of the likelihood 
of the Higgs search experiments normalized to its value in the case of pure 
background, the so called likelihood ratio IZ, to further constraint the p.d.f. 
for the Higgs mass obtained employing only the precision physics data, 
f{mH I ind.). 

In this paper we would like to recall the main features of our method and 
present an updated analysis. With respect to Ref.||T| we improve our anal- 
ysis in several aspects: i) we use the exact TZ for searches up to y/s = 196 
GeV as provided by the LEP Higgs Working Group ii) We include as 
observables most sensitive to the Higgs boson mass not only the effective 
mixing parameter, s\v? O^^J*^ = Sgjj, and the W boson mass, Mw, but also 
the Z° leptonic partial width, F^. iii) We take into account small non lin- 
earity effects in the theoretical formulae, iv) We use the most recent results 
on the various precision observables. 

2 Role of the likelihood ratio in reporting results 
of searches 

We begin by discussing the role of the likelihood ratio in constraining the 
Higgs mass. Let us assume that through the information coming from pre- 
cision measurements we have obtained f{mH \ ind.). A natural question to 
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ask is then how this p.d.f. should be modified in order to take into account 
the knowledge that the Higgs boson has not been observed at LEP for cen- 
ter of mass energies up to the highest available. To answer this question 
let us discuss first an ideal case. We consider a search for Higgs production 
in association with a particle of negligible width in an experimental situ- 
ation of "infinite" luminosity, perfect efficiency and no background whose 
outcome was no candidate. In this situation we are sure that all mass values 
below a sharp kinematical limit Mk are excluded. This implies that: a) the 
p.d.f. for Mh must vanish below Mk; b) above Mk the relative probabili- 
ties cannot change, because there is no sensitivity in this region, and then 
the experimental results cannot give information over there. For example, if 
Mk is 110 GeV, then /(200 GeV)//(120 GeV) must remain constant before 
and after the new piece of information is included. In this ideal case we have 
then 

r mH < Mk 

f{mH\dir.&,ind.) = l ^'T'^IH]. m„ > Mk , ^ 

where the integral at denominator is just a normalization coefficient. 

More formally, this result can be obtained making explicit use of the 
Bayes' theorem. Applied to our problem, the theorem can be expressed as 
follows (apart from a normalization constant) : 

f{mH I dir. & ind.) oc f{dir. \ mn) ■ f{mH \ ind.) , (2) 

where f{dir \ mn) is the so called likelihood. In the idealized example we are 
considering now, f{dir \ m//) can be expressed in terms of the probability of 
observing zero candidates in an experiment sensitive up to a Mk mass for 
a given value itih, or 

f{dir. I niH) = /("zero cand." | niH) = | J mH>MK 

In fact, we would expect an "infinite" number of events if M^ were below 
the kinematical limit. Therefore the probability of observing nothing should 
be zero. Instead, for Mfj above Mk, the condition of vanishing production 
cross section and no background can only yield no candidates. 

Consider now a real life situation. In this case the transition between 
Higgs mass values which are impossible to those which are possible is not 
so sharp. In fact because of physical reasons (such as threshold effects and 
background) and experimental reasons (such as luminosity and efficiency) 



3 



we cannot be really sure about excluding values close to the kinematical 
limit, nevertheless the ones very far from Mx are ruled out. Furthermore, 
the kinematical limit is in general not sharp. In the case of Higgs produc- 
tion at LEP the dominant mode is the Bjorken process e~^e~ — > H + Z° . 
Indeed, this reaction does not have a sharp kinematical limit at y/s — Mz 
(minus a negligible kinetic energy), due to the large total width of the Z° . 
The effective kinematical limit [Mx^^f) depends on the available integrated 
luminosity and could reach up to the order of ^ y/s — Mz + 0(10 GeV) for 
very high luminosity. Thus, in a real life situation we expect the ideal step 
function likelihood of Eq. (^ to be replaced by a smooth curve which goes 
to zero for low masses. Concerning, instead, the region of no experimental 
sensitivity, Mh ^ Mx^ff, the likelihood is expected to go to a value inde- 
pendent on the Higgs mass that however is different from that of the ideal 
case, i.e. 1, because of the presence of the background. 

In order to combine the various pieces of information easily it is con- 
venient to replace the likelihood by a function that goes to 1 where the 
experimental sensitivity is lost Q. Because constant factors do not play any 
role in the Bayes' theorem this can be achieved by dividing the likelihood 
by its value calculated for very large Higgs mass values where no signal is 
expected, i.e. the case of pure background. This likelihood ratio, TZ, can be 
seen as the counterpart, in the case of a real experiment, of the step function 
of Eq. (^). Therefore, the Higgs mass p.d.f. that takes into account both 
direct search and precision measurement results can be written as 

/(m„ I <i.r. & »d) = .o.i^'"'"'^"'" . (4) 

jg Tc{mH) J [rriH \ md.) Amu 

In Eq. (Q) 7^, namely the information from the direct searches, acts as a 
shape distortion function of /{mn \ ind.). As long as lZ{mH) is 1, the shape 
(and therefore the relative probabilities in that region) remains unchanged, 
while TZ^mn) indicates regions where the p.d.f. should vanish. A 
conventional limit can be derived by the TZ function alone in the transition 
region between the region of firm exclusion (TZ — > 0) and the region of 
insensitivity {TZ 1). However this limit can only have the meaning of 
a 'sensitivity bound' |^], and cannot be a probabilistic limit which tells us 
how much we are confident that the Higgs mass is above a certain value. To 
express consistently our confidence we need to pass necessarely through (Q) . 

One should notice that TZ^mn) can also assume values larger than 1 for 
Higgs mass values below the kinematical limit. This situation corresponds to 
a number of observed candidate events larger than the expected background. 
In this case the role played by TZ{mH) is to stretch /{mn \ ind.) below the 
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effective kinematical limit and this might even prompt a claim for a discovery 
if TZ becomes sufficiently large for the probability of Mh in that region to 
get very close to 1. 

3 Higgs mass inference from precision measure- 
ments 

We are going to construct /{rriH \ ind.) employing the three observables, 
Sgjj, Mw and T^. These quantities are the most sensitive to the Higgs mass 
and also very accurate measured. The most convenient way to approach the 
problem is to make use of the simple parameterization proposed in Ref.|^] 
and updated in Ref.|p where Sgjj, M]y and are written as functions of 
Mh, Mt, as and the hadronic contribution to the running of the electro- 
magnetic coupling: 

Se// = (Se//)o + Ci^i + C2^2 - C3^3 + C4^4 , (5) 

Mw = - diAi - d^Al - d2A2 + dsAs - d4A4 , (6) 

= r?- (71^1 -55^? -52 ^2 + 53 ^3 -54 ^4. (7) 

In the above equations Ai = ln(MH/100GeV), A2 = [(Aa)/,/0.0280 - 1], 
^3 = [(Aft/175 GeV)2 - 1] and A^ = [{as{Mz)/0.n8 - 1], where Mt is the 
top quark mass, as{Mz) is the strong coupling constant and {Aa)h is the 
five-flavor hadronic contribution to the QED vacuum polarization at q"^ = 

. (•Se//)o5 ^^^^ '^e excellent approximation) the theoretical 

results obtained at the reference point (Aa)/j = 0.0280, Mj = 175 GeV, 
and as{Mz) = 0.118 while the values of the coefficients Cj, di and gi are 
reported in Tables 3-5 of Ref.[^ for three different renormalization scheme. 
Formulae (|5[-^) are very accurate for 75 <^ Mh ^ 350 GeV with the other 
parameters in the ranges 170 <, Mt ^ 181 GeV, 0.0273 <, (Aa)^ <, 0.0287, 
0.113 ^ as{Mz) ^ 0.123. In this case they reproduce the exact results of 
the calculations of Refs.(||7|,[^,@) with maximal errors of Ss^j-jr ~ 1 x 10^^, 
SMw <, 1 MeV and 5Ti <^ 3 KeV, which are all very much below the 
experimental accuracy. Outside the above range, the deviations increase 
but remain very small for larger Higgs mass, reaching about 3 x 10~^, 3 
MeV, and 4 KeV at Mh = 600 GeV for sl^j, Mw, F^, respectively. 

Formula (^) can be seen as providing an indirect measurement of Ai = 
[sj.ff — (Se//)o — C2A2 + C3A3 — C/iA/i\/ci, while (|6|) and (^) of the quantities 
Y = M^-Mw- d2A2 + ^3^3 - ^4^4 and Z = F? - F^ - 52^2 + 53^3 - 
g^A^, respectively, all three variables being described by Gaussian p.d.f.'s. 



5 



In Ref.[|^ we considered only the relations (^) and (^) and used the fact 
that the non linearity effect given by the coefficient is small (in the MS 
scheme di = 5.79 • 10~^, d^ = 8.0 • 10~^) to linearize in Ai in order to 
directly obtain a determination also from M\y. Here, instead, we take 
into account exactly this non linearity effect, in view also of the fact that 
the uncertainty on Mw shrank. In Fig. |^ we plot the p.d.f. of Ai using 
formulae and the input quantities specified in the next section. The 

inference from s'^jf, M\y, and Ti, separately, and from their combination is 
shown. The quadratic expression in @ and (0) for Mw and Ti, respectively, 
gives rise to an unphysical peak on the left side of the plot for values of the 
Higgs mass where formulae (^-0) are not valid. Neverthless, as shown in the 
bottom plot, when the total combination is taken these unphysical peaks 
disappear and the inference obtained is concentrated in the physical region. 

The Ai, Y and Z determination are clearly correlated, therefore one has 
to built a covariance matrix. This can be easily done because formulae (^-@) 
are linear in the common terms X = {Mf, as{Mz), {Aa)h}- The likelihood 
of our indirect measurements Q = {Ai, Y, Z} is then a three dimensional 
correlated normal with covariance matrix 

--.^Ell-t-^^™ (B) 

or 

/(^l In(mH)) oce-^'/2 (g) 

where = A^V~^A with = {ai - ln(mH/100), y - di ln(mH/100) - 
d5ln2(m///100), z - 51 ln(m///100) - g^ln^imn /WO)}. Using Bayes' theo- 
rem the hkelihood @ can be turned into a p.d.f. through the choice of a 
prior. The natural choice is a uniform prior in ln(m//), since it is well under- 
stood that radiative corrections measure this quantity (for this reason the 
likelihood has been expressed in terms of ln(m//)). Moreover, this choice 
recovers the result of Ref. which was obtained as uncertainty propaga- 
tion without explicit use of a prior on the Higgs mass. The uniform prior in 
In(mii-) implies that /(ln(m/f/100) | ind.) is just the normalized likelihood 
(P). Using standard probability calculus we can express our results as a 
p.d.f. oiMn: 

f(raH I ind) = _^ -(yV2) a " ^^^^ 

The theoretical coefficients, Cj, dj, gi, depends on the renormalization 
scheme in which the relevant calculations are done and their numerical 
spread is usually taken as an indication of the theory uncertainty of the 
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Figure 1: Probability density function of Al = ln(Mj7/100 GcV) obtained 
from Sgjj, Mw, and their combination for different values of {Aa)h (see 
comment on the text about low mass peaks). 
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calculations. A way to take into account this uncertainty is to consider dif- 
ferent inferences, each conditioned by a given set of parameters, labelled by 
Ri. For each renormalization scheme Ri we can construct a f{mH \ ind., Ri) 
and obtain a p.d.f. "integrated" over the possible schemes, through 

/(m// I ind.) = ^ /(m// | ind.,Ri) ■ f{Ri) , (11) 

i 

where f{Ri) is the probability assigned to each scheme {f{Ri) = 1/3 Vz). 
The calculation of expectation value and variance is then straightforward or 

E[Mh] = ^Y.^[MH\Ri] (12) 

i 

ct\Mh) = ^J2^^(^^^i\R^) + lY.^^i^H\R^]-E^[MH] (13) 

i i 
i 

where aE indicates the standard deviation calculated from the dispersion of 
the expected values. In Ref. [l[ we have shown that results almost identical 
are obtained if one employs for the coefficients Cj, di an average value and 
a standard deviation evaluated from the dispersion of the values obtained 
from the various renormalization schemes. 



4 Results 

The experimental inputs we use to construct /(m/^ | ind.) are s^jj = 
0.23151 ± 0.00017, Mw = 80.394 ± 0.042 GeV, = 83.96 ± 0.09 MeV, 
Mt = 174.3 ± 5.1GeV, as{Mz) = 0.119 ± 0.003. Concerning {Aa)h in the 
recent years there has been a lot of activity on this subject with several 
evaluations. They can be classified as of two types: i) the most phenomeno- 
logical analyses, that rely on the use of all the available experimental data 
on the hadron production in e+ annihilation and on perturbative QCD 
(pQCD) for the high energy tail {E > 40 GeV) of the dispersion integral. 
The reference value in this approach is (Aa)f'^ = 0.02804 ± 0.00065 @. 



ii) The so called "theory driven" analyses |pT| , that differ from the pre- 
vious most phenomenological ones mainly by the use of pQCD down to 
energies of the order of a few GeV and by the treatment of old experi- 
mental data in regions where pQCD is not applicable. The combination 
of these two factors gives a result for (Aa)/j that differ from type i) one 
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by a drastically reduced uncertainty but, at the same time, a lower central 
value. The most stringent evaluation of these theory oriented analyses is 
(Aa)f^ = 0.02770 ± 0.00016 @, that we use as reference value for this 
kind of approach. At the moment there is no definite argument for choosing 
one or other of the two approaches. The results are absolutely compatible 
to each other. However, the numerical difference between central values 
and uncertainties is such that it prevents an easy estimation of the effect 
of choosing one value instead of the other. For these reasons we decided to 
present our results for the values of {Aa)h given by (Aa)^'^ and (Aa)^^ 
separately. 

The values of the TZ function that enters in Eq. has been provided by 
the LEP Higgs Working Group |2| : they take into account the Higgs searches 
by all four LEP collaborations for center of mass energy up to ^/s = 196 
GeV. 

Table || summarizes the result of our analysis in terms of various con- 
venient parameters of the distribution. We present the two cases {Aa)h = 
{Aa)^'^ and {Aa)h = {Aa)f^^ and report all values in TeV to reduce the 
number of digits to the significant ones. The shape of the p.d.f. with and 
without the inclusion of the direct search information is presented in Fig. ^. 
From this figure one can notice that, in the case of f{ind. \ mn), the use of a 
higher central value for {Aa)h (i.e. (Aa)^"^) tends to concentrate more the 
probability towards smaller values of Mh- As a consequence the analysis 
based on (Aa)f gives results for the standard deviation of the Mh p.d.f. 
and 95% probability upper limit, Mfj, very close to those obtained using 
(Aa) regardless the fact that the uncertainty on (Aa)^"' is approximately 
4 times larger than that of {Aa)f^^ . As expected, the inclusion of the direct 
search information in the Higgs mass probability analysis drifts the p.d.f. to- 
wards higher values of Mh, changing its shape such that the probability of 
Mh values below 110 GeV drops to « 9%. The various parameters of the 
distribution (expected value, standard deviation, mode (M^) and median 
(Mfj)) are not very sensitive to the values of the hadronic contribution to 
the vacuum polarization. Also in both cases, ~ 80 % of the probability is 
concentrated in the region Mh < 0.20 TeV. Instead the choice of {Aa)h 
affects the tail of the distribution with (Aa)^"^ producing a much longer 
one. 
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{Aa)h = 0.02804(65) 



{Aa)h = 0.02770(16) 



E[MH]/TeV 
cj(MH)/TeV 
Mn/TeV 
M^°/TeV 

P{Mh < 0.11 TeV) 
P{Mh < 0.13 TeV) 
P{Mh < 0.20 TeV) 

M^/TeV; P{Mh < M^f) 
Mfi^/TeV; pImh < Mf ) 

jMi/TeV; P{Mh < Mi) 
1 Ma/TeV; P{Mh > Ma) 



! 0.95 
! 0.99 

0.16 
0.16 



{ind. 



0.10 
0.07 
0.06 
0.09 

67% 
76% 
92% 

0.23 
0.33 





8.9% 
34% 
79% 

0.29 
0.41 

/0.12 
1 0.22 



{ind. 



0.12 
0.06 
0.09 
0.11 

51% 
65% 
91% 

0.23 
0.31 




ind. 

+ 

y dir. 

0.16 
0.05 

0.11 
0.14 

8.5% 

35% 
83% 

0.26 
0.34 

f 0.12 
I 0.20 



Table 1: Summary of the direct plus indirect information. 



5 Conclusions 

The likelihood ratio IZ \s a, form to report the experimental results that 
allows an easy combination of the information from Higgs search experiments 
with that coming from precision measurements and accurate calculations in 
order to constraint jointly the Higgs mass. The TZ function is very convenient 
for comparing and combining the various informations and it has also an 
intuitive interpretation because of its limit to the step function of the ideal 
case. It can be also seen as the p.d.f. for the Higgs mass in the case of 
complete lack of other information, i.e. when one assumes f{m.H \ ind.) = 1. 
Although in this case the normalization integral in Eq. (^) is mathematically 
"infinite" , still the relative probabilities of different intervals of mass regions 
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Figure 2: Probability distribution functions using only indirect information (solid 
line) and employing also the experimental results from direct searches (dashed 
one): a) {Aa)h = 0.02804(65); b) (Aa);, = 0.02770(16). 

are perfectly well defined. However, it is obvious that it is not possible to 
evaluate from the TZ function alone a probabilistic lower limit. This can 
only be done when TZ is combined with fimn \ ind.) obtained from precision 
measurements, which has the important role of making large values of Mh 
impossible, thus making the final p.d.f. normalizable. 

The analysis we have performed clearly shows that a heavy Higgs sce- 
nario is highly disfavored, the data preferring a \ogiQ{MH /GeV) k, 0{2) . 
Note that our results are derived under the assumption of the validity of the 
SM and rely on the input experimental and theoretical quantities stated in 
the text. All these assumptions seem to us very reasonable. In particular, 
we don't consider strong evidence against the SM the fact that about one 
half of f{mH I ind.) is eaten up by the LEP direct search. 



We wish to thank the LEP Higgs Working Group for presenting the likeli- 
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hood ratio values of the Higgs searches and P. Igo-Kemenes and G. Ganis 
for useful communications. 
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